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Plasma  Instability  and  Terahertz  Oscillations  in 
Resonant-Tunneling  Transistors 

Victor  Ryzhii  and  Michael  Shur 


Abstract — The  self-excitation  of  plasma  oscillations  in 
a resonant-tunneling  transistor  structure  is  studied.  It 
is  demonstrated  that  it  can  give  rise  to  the  generation  of 
transient  current  at  frequencies  in  the  terahertz  range. 
The  amplitudes  and  the  frequencies  of  different  modes 
of  the  generated  oscillations  are  voltage-tuned. 

resonance. 

I.  Introduction 

In  this  paper,  we  study  the  instability  of  standing 
plasma  waves  in  a resonant-tunneling  (RT)  transistor 
considering  the  nonlinear  stage  of  their  excitation  and 
the  generation  of  terahertz  radiation  at  fundamental 
plasma  frequency  and  harmonics.  The  structure  un- 
der consideration  consists  of  an  undoped  quantum  well 
sandwiched  by  a double-barrier  RT  structure  from  one 
side  and  thick  barrier  layer  from  another  side.  The 
QW  is  supplied  with  the  side  contacts.  This  structure 
is  sandwiched  between  narrow  gap  doped  layers  serv- 
ing as  the  collector  contact  and  the  gate,  respectively. 
The  QW  plays  a role  of  the  emitter  channel  in  which  a 
two-dimensional  (2D)  electron  gas  is  induced  by  bias- 
ing voltages  applied  between  the  emitter  channel  side 
contacts  and  both  the  collector  contact  and  the  gate. 
Similar  structures  were  proposed  and  realized  as  work- 
ing devices  earlier  [1],  [2].  However,  contrary  to  the 
normal  operation  of  these  devices,  we  assume  that  the 
polarity  and  the  value  of  the  emitter-collector  voltage 
( Vc ) are  chosen  to  provide  the  RT  current  from  the 
emitter  channel.  The  emitter-gate  voltage  (V?)  deter- 
mines the  electron  concentration  in  the  QW  emitter 
channel  and,  in  turn,  the  plasma  wave  velocity  and  fre- 
quencies. [3],  [4],  [5],  [6]  The  conduction  band  profile 
of  the  RT  structure  under  applied  voltages  is  shown  in 
Fig.  1.  Alternatively,  we  can  use  a structure  similar  to 
a standard  high-electron-mobility  transistor  but  with 
a double  barrier  RT  structure  placed  into  the  wide  gap 
region,  which  separates  the  channel  from  the  gate. 

II.  Model 

We  assume  that  the  2D  electron  gas  in  the  QW  emit- 
ter channel  is  degenerate.  The  spatio-temoral  varia- 
tions of  electron  sheet  concentration  £ = £(£,#),  aver- 
age lateral  velocity  in  this  channel  channel  u — #), 

and  potential  of  the  latter  <p  = p(t,  rr)  (with  respect 
to  the  QW  emitter  channel  contacts)  are  described  by 
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the  following  set  of  equations: 


dT,  dYlu 
dt  dx 
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Here  e and  m are  the  electron  charge  and  mass,  r"1 
is  the  product  of  the  try-to-escape  frequency  and  the 
maximum  transmission,  v is  the  electron  collision  fre- 
quency, Wc  and  Wg  are  the  thicknesses  of  the  RT 
structure  and  the  emitter-gate  barrier,  £*  is  the  in- 
duced equilibrium  concentration  of  electrons  deter- 
mined by  their  threshold  voltage  and  the  bias  volt- 
ages, £rt  = £rt  — 2eaiP  anc^  ^ are  energy  , and 
the  width  of  RT  resonance,  respectively,  and  a ~ 1 is 
the  geometrical  factor.  The  axis  x is  directed  in  the 
emitter  channel  plane.  The  RT  resonance  form-factor 
is  approximated  as  A(z)  — (1  + 22)-1.  Equation  (3) 
is  a consequence  of  the  Poisson  equation  in  the  so- 
called  weak  non-locality  approximation.  This  equation 
is  more  general  than  that  often  used  and  which  corre- 
sponds to  the  “gradual  channel”  approximation  [7]. 

We  shall  further  assume  that  the  span  of  the  QW 
emitter  potential  variations  is  not  too  large  that  non- 
linearities  associated  with  the  left-hand-side  terms  of 


Fig.  1.  Conduction  band  profile  of  the  RT  structure. 
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Eqs.  (1)  and  (2),  i.e.,  hydrodynamic  nonlinearities,  are 
weak.  When  the  width  of  the  tunneling  resonance 
is  small,  so  that  this  resonance  is  sufficiently  sharp, 
the  hydrodynamic  nonlinearities  can  be  neglected  in 
comparison  to  the  nonlinearity  caused  by  the  volt- 
age dependence  of  the  RT  current  in  the  right-hand 
side  of  Eq.  (1).  Let  us  chose  such  a value  of  volt- 
age Vc  to  provide  Srt  = T/y/3.  This  voltage  corre- 
sponds to  negative  differential  conductivity  of  the  RT 
structure  and  the  maximum  of  its  modulus.  Under 
these  assumptions,  the  right-hand  side  of  Eq.  (1)  can 

be  expanded  in  powers  of  <p  and  presented  as  -[jo  — 

(crtP  + iRT'J?)],  where  jo  is  the  steady-state  compo- 
nent of  the  current  density,  cfrt  — ~(\/3eajo/4r),  and 
7 rT  = (3\/3/64)(ea/r)3<;o  is  the  parameter  of  nonlin- 
earity of  the  RT  structure. 

The  boundary  conditions  for  the  potential  of  emitter 
channel  are  set  at  the  side  contacts: 

<p\x  = ±L  = °>  (4) 

where  2 L is  the  channel  length  (see  Fig.  1).  In  this  case, 
searching  for  the  emitter  channel  potential  in  the  form 
<p(t,x)  = + cos  knx,  where  tpo(x)  is  the 

n 

steady-state  component  associated  with  the  steady- 
state  current  along  the  channel,  kn  = 7rn/2L  and  n 
is  the  mode  index  ( n — 1,2,3,...),  When  the  poten- 
tial potential  drop  along  the  channel  is  small,  i.e., 
ea\^po(x)\  T (see  the  estimate  below),  the  nonuni- 
formity of  the  steady-state  components  of  potential, 
electron  concentration  and  current  density  can  be  ne- 
glected (<£>o  ^ 0,  So  — = const , and  jo  ~ const), 

and  Eqs.  (1)  - (3)  give  rise  to  the  following  equation: 

[(^+£)("-+£)+fcH^ 

-~(v  + Jl)  T (5) 

Here  vn  = v0/bn,  sn  = s0/y/K.,  and  $£n"n'"  = 
(47r/^e)7/^rW<5",n"n",  /bn,  where  vq  = (An  / &)<7  rtW  — 
-(A/37r/ae)(eaW/r)io,  «o  = (47re2S0W//aem)1/2,  is 
the  plasma  wave  velocity  in  the  limit  of  small  kn, 
bn  = 1 + (k2nWcWg/ 3),  W = WcWg/(Wc  + Wg), 
and  6™  71  n"'  are  the  coefficients  determined  by  the 
overlapping  integrals  for  different  interacting  plasma 
modes.  In  particular,  = 1 and  0\11  = 6\,  where 
a __  ( ea\3  Wjo 

1 *“  v 16a e )\T  ) 1 + (*2WcWg/12L2Y 

III.  Instability 

One  can  show  that  the  state  with  E = So  and  ip  — 

0 can  be  unstable  with  respect  to  the  excitation  of 
plasma  oscillations.  Indeed,  in  the  most  interesting 
case  when  k^s^  + v vn  > 0,  assuming  that  is  small 
and  neglecting  the  nonlinear  term  in  the  right-hand 


side  of  Eq.  (5),  one  can  arrive  at  the  following  formulas 
for  the  frequency  of  plasma  oscillations  (u>n)  and  their 
growth  rate  (7n): 


OJn  ~ 


Ksl- 


{y  - vn) 


1/2 


In  = 


{y  + vn) 


(6) 


One  can  see  from  Eqs.  (6)  that  when  the  differen- 
tial conductivity  of  the  RT  structure  is  negative  and, 
hence,  uq  < 0,  the  uniform  state  under  consideration 
can  be  unstable  (7 n > 0)  with  the  excitation  of  differ- 
ent plasma  modes  having  the  frequencies  un.  In  this 
case,  the  criterion  of  instability  has  the  form  v+vn  < 0. 
The  latter  can  be  presented  as 

i/  3 L2  ' V ' 

As  can  be  seen  from  this  inequality,  the  instability  oc- 
curs if  the  differential  conductivity  is  negative  and  its 
modulus  is  sufficiently  large.  Inequality  (7)  yields  the 
instability  criterion  different  from  that  obtained  previ- 
ously. [8]  The  difference  is  due  to  the  second  term  in 
the  right-hand  side  of  inequality  (7)  which  reflects  the 
effect  of  nonlocality.  Such  an  effect  imposes  a restric- 
tion on  the  excited  plasma  modes  if  the  ratio  \vo\/v  \§ 
fixed.  Thus,  in  the  situation  under  consideration  the 
spectrum  of  the  excited  plasma  oscillations  is  limited 
to  relatively  small  mode  indices  n.  An  additional  sup- 
pression of  the  instability  of  the  modes  with  large  in- 
dexes can  be  due  to  the  nonideality  of  the  2D  electron 
gas  in  the  channel  and  due  to  the  delay  of  electrons 
associated  with  the  RT  tunneling. 


Fig.  2.  Frequencies  of  first  and  second  modes  versus  channel 
length  for  n = 5 X 104  cm2/Vs. 


Figure  2 shows  the  frequencies  of  the  different 
plasma  oscillation  modes  as  function  of  the  channel 
length  for  different  electron  concentrations,  for  the 
structure  with  Wc  = 2 x 1CP6  cm,  The  structural 
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Fig.  3.  Growth  rate  of  different  plasma  modes  versus  chan- 
nel length  at  different  electron  concentrations  for  fi  = 5 X 
104  cm2/Vs. 


parameters  used  for  this  and  subsequent  figures  are 
as  follows:  Wg  = 4 x 10"6  cm,  T — 7.5  meV,  and 
re  = 17.5  ps,  One  can  see  that  at  reasonable  device 
sizes  and  electron  concentrations,  the  plasma  frequen- 
cies fall  into  the  terahertz  range.  Figure  3 shows  the 
dependencies  of  the  plasma  oscillations  growth  rate 
versus  the  channel  length,  As  can  be  seen  from  Fig.  3, 
the  increase  in  the  electron  concentration  leads  to  a 
significant  rise  of  the  growth  rate  of  plasma  oscilla- 
tions and  the  marked  increase  of  the  range  where  such 
a growth  occurs.  The  use  of  the  channel  with  higher 
electron  mobility  favors  the  instability  conditions. 


IV.  Nonlinear  oscillations 

We  shall  search  for  the  solution  of  Eq.  (5)  in  the 
form  ipn(t)  = An(t)  cos[<jjnt  — an(£)],  where  amplitudes 
An{t)  and  phases  an(t)  are  relatively  slowly  varying 
functions  of  time.  Further  we  shall  assume  that  only 
the  fundamental  mode  (n  = 1)  of  the  plasma  oscil- 
lations has  a positive  increment.  Substituting  i/jn  into 
Eq.  (5)  and  averaging  over  the  period  of  the  plasma  os- 
cillation, near  the  instability  threshold,  where  \v  + i/±\ 
is  small,  one  can  obtain  the  following  equation  for  the 
amplitude  of  fundamental  mode  (n  = 1)  of  plasma  os- 
cillations (compare  with  [9]): 


dA i 
dt 


o A1 


(8) 


da  1 _ {v_  + </i)2  3 {v\  - v)9i  2 

dt  4u>i  16  wi  11 


(9) 


Deriving  Eqs.  (8)  and  (9),  we  neglected  the  nonlinear 
effect  of  higher  modes  (n  > 1)  because  such  modes 
are  assumed  to  be  stable  in  the  linear  approximation. 
However,  the  excitation  of  the  fundamental  mode  can 


result  in  the  forced  excitation  of  the  stable  modes. 
Nevertheless,  when  |iq|  barely  exceeds  v , one  can  as- 
sume that  the  amplitudes  of  higher  modes  are  rather 
small. 

As  follows  from  Eq.  (9),  the  temporal  evolution  of 
the  amplitude  of  the  oscillations  is  given  by 

/ ~ A0  exp (-§(»  + yQt) 

“ 1 [1  + (A* /A?)2  exp(— (1/  + i/O*)] 1/2  ’ 

where  A?  is  the  initial  amplitude  (its  fluctuational 
value)  and  AJ°  is  the  final  steady-state  amplitude: 
AJ°  = [-4(v  + v^/SOi]1/2.  This  state  is  stable,  be- 
cause for  its  small  perturbations  Eq.  (9)  yields  = 
v + vx  < 0.  Near  the  stable  state,  the  oscillation 
phase  slowly  varies.  This  corresponds  to  the  fol- 
lowing nonlinear  correction  of  the  plasma  frequency: 
Auq/ui  ~ v{v  + Ui)/cjx. 

Using  the  obtained  amplitude  of  oscillations  of  the 
potential  in  the  QW  emitter  channel  one  can  arrive 
at  the  following  expressions  for  the  amplitudes  of  the 
first  ( 8J1  oc  Ai/F)  and  third  (8J3  oc  A3/F3)  harmon- 
ics of  the  electron  current  from  the  QW  emitter  chan- 
nel to  the  collector  contact  caused  by  the  fundamental 
plasma  harmonic: 


8J1 

Jo 


(>-f 


Jth\1/2  6J3 


OC 


where  J*‘ = ( H)  (w)  (£) 


and  Jq  = 2Ljo  is 


the  dc  emitter-collector  current.  The  amplitudes  of  ac 
electron  current  can  also  be  expressed  via  the  electron 
concentration  in  the  QW  emitter  channel:  5Jx/Jo  oc 
[1  - (E^/Sq)]1/2  and  SJ3/J0  oc  [1  - (Et/l/E0)]3/2, 

where  — —7= — 0 *T.  Going  back  to  Eq.  (5),  one 

V37 r ezaW 

can  find  that  the  coupling  between  the  fundamental 
(unstable)  and  higher  (stable)  modes  should  not  result 
in  a marked  excitation  of  the  latters  if  the  amplitude  of 
fundamental  mode  is  not  too  large,  i.e.,  near  its  exci- 
tation threshold:  (Jo  - Jth)  < Jth  or  (E0  - %th)  < E ih 
Figure  4 shows  the  relative  amplitude  SJi/Jo  of  the 
established  current  oscillations  at  the  fundamental  fre- 
quency as  a function  of  the  electron  concentration  in 
the  QW  emitter  channel  for  different  electron  mobil- 
ities in  the  latter.  Due  to  the  dependence  of  the  in- 
duced electron  concentration  in  the  QW  emitter  chan- 
nel on  the  bias  voltages,  particularly,  the  emitter-gate 
voltage,  E0  a Vg  — Vth,  (where  Vth  is  the  threshold 
voltage),  the  amplitude  of  ac  electron  current  can  be 
effectively  controlled  by  this  voltage.  The  emitter-gate 
voltage  can  also  tune  the  frequency  of  the  excited  cur- 
rent oscillations. 

Using  Eqs.  (1)  - (3),  one  can  estimate  the  potential 
drop  along  the  channel  as  max  e\(po\  ~ JoiymL  / e&o) . 
Accounting  for  this  estimate,  one  can  find  that  con- 
ditions max  ea\ipo(x)\  <C  F (weakness  of  the  steady- 
state  potential  nonuniformity)  and  Jo  > Jth  (see  re- 
lationships (11))  do  not  contradict  to  each  other  if 
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Fig.  4.  Amplitude  of  oscillations  versus  electron  concentration 
in  the  channel  in  structures  with  different  electron  mobilities 
and  2 L — 0.3  /mi . 


v < s/L  ~ lj\.  This  inequality  corresponds  to  suf- 
ficiently high  electron  mobility  and  concentration  in 
the  channel  and  sufficiently  small  length  of  the  latter. 
It  is  typical  for  different  devices  utilizing  the  plasma 
waves  (see,  for  example,  [4],  [5],  [6]).  Consequently, 
the  electron  drift  should  not  lead  to  a marked  change  in 
the  instability  conditions.  However,  at  a large  poten- 
tial drop  along  the  channel,  one  may  expect  a marked 
change  in  the  RT  transistor  behavior  due  to  a combi- 
nation of  the  considered  mechanism  and  the  instability 
mechanisms  studied  previously. 

V.  Conclusions 

Using  an  analytical  model,  we  studied  the  self- 
excitation of  terahertz  current  oscillations  and  their 
nonlinear  stage  in  a RT  transistor  structure  associated 
with  the  plasma  instability  in  the  QW  emitter  chan- 
nel. The  criteria  of  the  instability  were  found  and  the 
oscillations  frequency  and  amplitude  were  calculated 
as  functions  of  the  structural  parameters  and  biasing 
voltage. 
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